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ABSTRACT 

Local luminous infrared (IR) galaxies (LIRGs) have both high star formation rates (SFR) and a 
high AGN (Seyfert and AGN/starburst composite) incidence. Therefore, they are ideal candidates to 
explore the co-evolution of black hole (BH) growth and star formation (SF) activity, not necessarily 
associated with major mergers. Here, we use Spitzer/IRS spectroscopy of a complete volume-limited 
sample of local LIRGs (distances of < 78Mpc). We estimate typical BH masses of 3 x 10 7 M Q 
using [Ne in] 15.56 /imand optical [Oiii]A5007 gas velocity dispersions and literature stellar velocity 
dispersions. We find that in a large fraction of local LIRGs the current SFR is taking place not only 
in the inner nuclear ~ 1.5 kpc region, as estimated from the nuclear 11.3 fim PAH luminosities, but 
also in the host galaxy. We next use the ratios between the SFRs and BH accretion rates (BHAR) 
to study whether the SF activity and BH growth are contemporaneous in local LIRGs. On average, 
local LIRGs have SFR to BHAR ratios higher than those of optically selected Seyferts of similar AGN 
luminosities. However, the majority of the IR-bright galaxies in the RSA Seyfert sample behave like 
local LIRGs. Moreover, the AGN incidence tends to be higher in local LIRGs with the lowest SFRs. 
All this suggests that in local LIRGs there is a distinct IR-bright star forming phase taking place prior 
to the bulk of the current BH growth (i.e., AGN phase). The latter is reflected first as a composite 
and then as a Seyfert, and later as a non-LIRG optically identified Seyfert nucleus with moderate SF 
in its host galaxy. 

Subject headings: galaxies: nuclei — galaxies: Seyfert — infrared: galaxies 



1. INTRODUCTION 

One of the most fundamental relations in extragalac- 
tic astronomy is that, at least in the local universe, the 
masses of supermassive black holes (BH, with masses 
Mbh > 10 6 Af Q ) correlate with the stellar mass and 
velocity disp ersion of the bulges of t heir host galax- 
ies (see e.g. iMagorrian et al. I 119981: iGebhardt et al. I 
[2rMlMa"rconi fc Hunt Il2f)03t IHaring fc Rix 112004 This 
seems to imply that bulges and supermassive BH evolve 
together and regulate each other. 

Mergers of gas-rich galaxies are efficient in both pro- 
ducing elevated star formation rates (SFR) and trans- 
porting gas to the nuclear region to allow for BH growth. 
If sufficient matter becomes available very close to the nu- 
clear BH and is accreted, the nucleus of the galaxy will 
shine wi th enormous power as an active galactic nucleus 
(AGN, iLvnden-Bell I Il969i) . In parallel, and probably 
before the fully developed AGN phase, the merger will 
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trigger a high rate of star format ion leading to a lum i- 
nous infrared (IR)-active phase ([Sanders et al. l Fl988). 
either as a luminous IR galaxy (LIRG) or ultraluminous 
IR galaxy (ULIRG). These are defined as having an IR 
8-1000 ^m lum inosity L m > 10 11 or L m > 10 12 L Q , 
respectively (see lSanders fc Mirabel 1119961 for a review). 

As an alternative to this sequence. iKormendv et aL~l 
(|2011l) proposed that BHs in bulgeless galaxies and in 
galaxies with pseudobulges grow through phases of low- 
level Seyfert-like activity. Such growth is believed to 
be driven stochastically by local processes (secular pro- 
cesses), and thus it would not ha ve a global impact o n 
the host galaxy structure (see also lHopkins et aLl feOOS). 
Moreover, there is now evidence of two fundamentally 
different modes of BH growth at work in early- and late- 
type galaxies in the local universe (jSchawinski et al~1 
I2010T ). The role of moderate luminosity AGN that reside 
in late-type galaxies and do not involve any recent major 
mergers is not clear and needs further investigation. 

LIRGs are powered by both AGN an d sta r 
formation activity (ISanders &: Mirabel I Il996h . 
lAlonso-Herrero et al. I (|2012t ) showed that while an 
AGN accompanies star formation activity in a large 
proportion of local LIRGs, in most cases the AGNs are 
not energetically important. This is because the AGNs 
hosted in local LIRGs have Seyfert-like luminosities 
and thus their bolometric contribution to the total 
IR luminosity is small (typically 5%). Therefore, the 
IR luminosities of LIRGs imply SFRs in the range 
11 - HOMgyr" 1 using the iKennicuttl (j!998[) prescrip- 
tion converted to a lKroupal (|201 2T ) Initial Mass Function 
(IMF). 
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TABLE 1 




Spitzer/IRS Calibration sources 


Name 


AOR 


PNG043.1+37.7 


4109056 


PNG011. 7-00.6 


4109568 


PNG025.8-17.9 


4110080 


PNG265.7+04.1 


4111104 


PNG206.4-40.5 


4111616 


PNG285.7-14.9 


4112128 


PNG054.1-12.1 


4112640 


PNG009.4-05.0 


4114176 


PNG342.1+10.8 


4115200 


PNG316.1+08.4 


4116224 


Calwav-22B-NGC7027 


15343104 


IRSS-SEPN-N7293-0001 


15752960 



It is reasonable to assume that the same gas that is 
used to form stars in the host galaxy can also be used 
to feed the AGN provided that there is a mechanism 
able to transport the gas to the inner region (on scales 
of less than 0.1 pc) of the g alaxy (see e.g. the review of 
I Alexander fc Hickox I \201$) . Indeed, it is now apparent 
that there is a relation between star formation activity on 
sub-kiloparsec and bulge scales and the BH accretion rate 
(BHAR), as shown by observations and numerical sim- 
ulations dHeckman et al. I l2004t iHqpkins fc Quataert I 
I2010t iDiamond-Stanic k Rieke N2012D . This relation is 
also predicted to be present, although with a lower 
significance, for the integ rated SFR of the galaxy 
([Hopkins fc Quataert 1 120101 ). Therefore, LIRGs show 
propitious conditions to study the co-evolution of star 
formation activity and BH grow th as they show both 
high integrated and nuclear SFRs ()Alonso-Herrero et al. I 
2006) and a hi gh occurrence of AGN and composite nu- 
clei ( 60 - 70%. lYuan et al. 1120101 lAlonso-Herrero et al. I 
120121 ). 

In this paper we study the co-evolution of the SFR 
activity and the BH growth in the complete volume- 
limited sam ple of local LIRGs of lAlonso-Herrero et al. I 
(l2lMl2l)12l). We use Spitzer Infrared Spectrograph (IRS, 
iHouck et a l. 2004) observations to estimate the nuclear 
(~ 1.5 kpc) SFRs and compare them with the integrated 
values from the IR luminosities. We look for spectro- 
scopically resolved [Ne in] 15.56 fi,m lines. Assuming that 
this line is produced in the narrow line region (NLR) of 
the AGN, we measur e its velocity dispersion and use the 
iDasvra et al. I ()201lD relations to obtain the masses of 
the BH hosted in local LIRGs. We complement these 
with observations of the optical line [O in] A5007 and lit- 
erature values of the stellar velocity dispersion er* to ob- 
tain further estimates of the BH masses in local LIRGs. 
All this information allows us to determine the ratios 
between the SFR (both, nuclear and integrated) and the 
BHAR in local LIRGs and compare them with those of 
optically selected Seyfert galaxies. This comparison lets 
us explore the role of the LIRG-phase in the growth of 
BH in the local universe. Throughout this paper we as- 
sume the following cosmology Ho = 70 km s _1 Mpc _1 , 
ft M = 0.3 and 0\ = 0.7. 

2. SAMPLE, OBSERVATIONS, AND DATA REDUCTION 

2.1. The sample of local LIRGs 
For this work we use the complete volume-limited sam- 
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[NeIII]15.56/!zm FWHM (/im) 

Fig. 1. — Distribution of the observed (not corrected for instru- 
mental resolution) FWHMs of the [Ne III] 15. 56 £tm line of local 
LIRGs compared with those of the calibration sources. The latter 
are assumed to provide the instrumental resolution at the wave- 
length of the line. The dot represents the average value of the 
FWHM of the calibration sources with the corresponding 1<t error. 
We marked those LIRG nuclei classified as AGN (Seyfert and/or 
[Nev] emitters) and composites. The rest are Hn-like or have 
no classification. The thick line histogram represents the observed 
FWHM of those galaxies deemed to have spectroscopically resolved 
lines. 



pie of local LIRGs that was presented and discussed in selection criterion. 



detail by lAlonso-Herrero et al. I (|2006l l2012t ). We drew 
the samp le from the IRAS Rev ised Bright Galaxy Sample 
(RBGS, iSanders et al~ll200l to include all the sources 
with log(i m /L ) > 11.05 and v hel = 2750 - 5300 km 
s -1 . The sample is composed of 45 IRAS systems, with 
8 containing multiple nuclei. Therefore the sample in- 
cludes 53 individual galaxies. For the assumed cosmol- 
ogy the distances are in the range ~ 40 — 78 Mpc, with a 
median value of 65 Mpc. The IR luminosities of the indi- 
vidual galaxie^] are in the range log(LiR,/L0) = 10.64 — 
11.67, with a median value of log(LiR,/Lo) = 11.12. All 
the releva nt information for the s ample can be found in 
table 1 of lAlonso-Herrero et al~1 (|2012ft . 

Combining a number of optical and mid-in frared (mid- 
IR) indicators. lAlonso-Herrero et al. I (|2012f ) derived an 
AGN detection rate of ^62% for this complete volume- 
limited sample of local LIRGs. The derived AGN bolo- 
metric luminosities, from the mid-IR spectral decom- 
position and/or X-ray observations, are in the range 
Lboi(AGN) = (0.4 - 50) x lO^ergs" 1 with a median 
of ~ 1.4 X lO^ergs" 1 dPereira-Santaella et al. 1 120111: 
lAlonso-Herrero et al. I I2012T ) . However, these AGN are 
overall only responsible for ~ 5% of the IR luminosity 
emitted by local LIRGs, although the AGN contribution 
varies from source to sourc e (see llmanishi et al~l 120101 : 
lAlonso-Herrero et al. 1 12 012). This is in good agreement 
with iPetric et al. I (|2011l ) for The Great Observat ories 
All-Sky LIRG Survey (GOALS; lATmus et al. ll2009h . 

2.2. Spitzer/IRS Observations 

2.2.1. Observations 

We retrieved Spitzer/IRS (jHouck et al. I I2004D spec- 
troscopy of the sample of LIRGs taken with the high 

9 Galaxies significantly below the selection luminosity are mem- 
bers of multiple systems with total luminosities that satisfy the 
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resolution (R ~ 600) short-high (SH) and long-high 
(LH) modules that cover the 9.9 — 19.6 /im and 18.7 — 
37.2 (an spectral ranges, respectively. Details on the pro- 
gramme IDs, the observations, and data re duction of the 
galax ies, including t hose from GOALS (|Armus et al~l 
2009) , are given in [Alonso-Herrero et al. I ()2012f ) and 
Pereira-Santaella et al. I (|2010allbf ). We note that for var- 
ious reasons two galaxies in this volume-limited sample 
were not observed with the Spitzer/IRS high spectral 
resolution modules. We extracted the nuclear spectra 
assuming a point source calibration. 

We also retrieved SH and LH spectroscopy from the 
Spitzer archive of 12 calibration sources and planetary 
nebulae to obtain an accurate estimate of the spectral 
resolution using fine structure lines (see Section I2.2.3[) . 
The names and corresponding Astronomical Observation 
Request (AOR) numbers are give in Table [1] We reduced 
the data as for the sample of local LIRGs. 

2.2.2. Line Fitting 

For the fine-structure lines we used Gaussians to mea- 
sure the line flux, equivalent width (EW) , and full width 
half maximum (FWHM) and a first order polynomial 
to fit the local continuum (see iPereira-Santaella et al. I 
[20T0al for more details). The fluxes of the 
[Ne II] 12.81 fim, [Ne ml 15.56 /an, and [Oiv]2 5 89 /im 
line s can be found in IPereira-Santaella et al. I (|2010bT > 
and I Alonso-Herrero et al. I (|2012t ). For the 11.3 /im poly- 
cyclic aromatic hydrocarbon (PAH) feature first we fit- 
ted the local continuum using a linear fit between 10.6 
and 11.8 /tin. After subtracting the continuum, we in- 
tegrated the flux in the wavelength range between 10.8 
and 11.6 /im. We list the fluxes, EW, and corresponding 
errors of the 11.3 /im PAH feature in Table [2] 

2.2.3. Spectrally resolved [Nelll]15.56 [im lines 

The velocity dispersion of the ionized gas in the NLR 
allows measuring the BH mass (Mbh) in an AGN supple- 
m enting reverberation m apping techniques. As discussed 
bv lGreene &: Ho 1 (|2005f) . the NLR is sufficiently compact 
to be illuminated by the AGN, while large enough to 
feel the gravitational potential of the bulge. Therefore, 
there is a relatively good correlation between the gas and 
the stel lar velocity dispersion arou nd local AGN. In the 
mid-IR IDasvra et ai~l (|2008l 1201 lft showed that the ve- 
locity dispersion of the fine structure lines [S iv] 10.51 /im, 
[Ne ill] 15. 56 /im, [Nev] 14.32 /im, and [Oiv]25.89 /im of 
AGN are well correlated with -Mbh- 

For our Spitzer/IRS spectra we focus on the 
[Ne in] 15.56 /xm line to look for spectrally resolved line 
profiles. In AGN the flux and luminosity of this fine- 
structure line are found to correlate well with those of the 
mid-I R [Nev] lines at 14.3 and 24.3 / xm (|Goriian et al. I 
120071 IPereira-Santaella et al. 1l2010bh . The high ioniza- 
tion potential of the [Nev] lines indicates that they are 
mostly excited by AGN and thus, the [Ne in] 15.56 /im 
emission is likely to be as well. When compared with 
other mid-IR fine structure lines with relatively high ion- 
ization potentials (e.g., [Siv]10.51 /im, [0iv]25.89 /im, 
and the [Nev] lines), the [Ne in] 15.56 /im line is a good 
compromise between ionization potential, critical den- 
sitj{_l intensity of the line, and a clean region of the 

10 p or re f erence the [Ne III] 15. 56 (im line has a slightly higher 



TABLE 2 

Fluxes and EW of the nuclear 11.3 /jm PAH feature 



Galaxy Flux EW 



NGC 23 


430.4 ± 8.5 


0.78 


± 





05 


MCG+12-02-001 


396.8 ± 7.9 


0.42 


± 





01 


NGC 633 


166.1 ± 4.0 


0.61 


± 





01 


ESO 297-G012 


128.4 ± 3.5 


0.60 


± 





05 


UGC 01845 


347.5 ± 5.2 


0.74 


± 





03 


UGC 02982 


379 ± 35 


0.84 


± 





15 


CGCG 468-002 NED01 


59.6 ± 5.6 


0.24 


± 





03 


CGCG 468-002 NED02 


78.4 ± 7.7 


0.46 


± 





07 


UGC 03351 


275.5 ± 8.3 


0.78 


± 





06 


NGC 2369 


357.1 ± 3.9 


0.58 


± 





02 


NGC 2388 


428 ± 14 


0.53 


± 





01 


MCG+02-20-003 


149.5 ± 4.1 


0.70 


± 





06 


NGC 3110 


299.0 ± 4.0 


0.80 


± 





03 


NGC 3256 


1547.2 ± 4.8 


0.556 


± 





004 


ESO 264-G057 


165 ± 13 


0.76 


± 





11 


IC 694 


446 ± 16 


0.38 


± 





02 


NGC 3690 


308 ± 25 


0.09 


± 





01 


ESO 320-G030 


450.3 ± 4.6 


0.78 


± 





03 


MCG-02-33-098 W 


99.0 ± 4.2 


0.27 


± 





02 


MCG-02-33-098 E 


102.3 ± 6.8 


0.62 


± 





12 


IC 860 


43.2 ± 2.2 


0.44 


± 





01 


MCG-03-34-064 


23.0 ± 3.2 


0.013 


± 





002 


NGC 5135 


376.9 ± 4.0 


0.46 


± 





01 


ESO 173-G015 


434 ± 31 


0.53 


± 





06 


IC 4280 


249.4 ± 7.0 


0.80 


± 





01 


UGC 08739 


146 ± 41 


0.84 


± 





38 


ESO 221-IG010 


259.9 ± 9.4 


0.63 


± 





05 


NGC 5653 


385 ± 10 


0.81 


± 





07 


NGC 5734 


402 ± 11 


0.80 


± 





06 


NGC 5743 


251 ± 14 


0.75 


± 





12 


IC 4518 E 


64.6 ± 1.4 


0.62 


± 





05 


IC 4518 W 


62.6 ± 5.9 


0.16 


± 





02 


Zw 049.057 


62.3 ± 4.5 


0.96 


± 





23 


NGC 5936 


240.0 ± 6.6 


0.67 


± 





05 


NGC 5990 


243 ± 56 


0.29 


± 





09 


NGC 6156 


141 ± 16 


0.22 


± 





03 


IRAS 17138-1017 


307.6 ± 4.3 


0.53 


± 





02 


IRAS 17578-0400 


177 ± 10 


0.76 


± 





09 


IC 4687 


478.2 ± 4.2 


0.70 


± 





02 


IC 4734 


195.1 ± 3.8 


0.59 


± 





03 


NGC 6701 


242.2 ± 4.2 


0.67 


± 





03 


MCG +04-48-002 


394 ± 33 


0.73 


± 





11 


NGC 7130 


230.5 ± 3.7 


0.36 


± 





01 


IC 5179 


391.9 ± 3.6 


0.68 


± 





02 


NGC 7469 


583.5 ± 7.7 


0.237 


± 





005 


NGC 7591 


142 ± 12 


0.36 


± 





06 


NGC 7679 


363 ± 16 


0.61 


± 





05 


NGC 7769 


85 ± 13 


0.57 


± 





21 


NGC 7770 


131 ± 18 


0.47 


± 





13 


NGC 7771 


359 ± 16 


0.65 


± 





07 



Notes. The fluxes are in units of 10 14 erg cm 2 s 1 and the 
EW are in units of fim. 



mid-IR spectrum for accurate measurements. 

Using the calibration sources (Section I2.2.ip we ob- 
tained twelve independent measurements of the unre- 
solved width of the [Ne in] 15.56 /im line and inferred an 
instrumental width of FWHM inst = 0.02618 ± 0.0007 /im. 
This is equivalent to a spectral resolution of R = 
A/FWHMi ns t ~ 595 ± 16 and an instrumental veloc- 
ity dispersion of <7i ns t = 215 ± 6kms~ at the wave- 

ionization potential than the optical [Om]A5007 emission line but 
the optical line has a slightly higher critical density (see figure 6 of 
IDasvra et "aTni20ill ) . 
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FlG. 2. — Observed profiles (solid lines and filled dots) of the spectrally resolved [Nelll]15.56 ^.m lines of local LIRGs that are classified 
as a Seyfert and/or are [Nev] emitters. The shaded area shows a Spitzer/lYiS SH un resolved profile represented as a Gaussian with 
FWHM= 0.02618 //m, as determined from Spitzer/IBS calibration sources (Section 
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FlG. 3. — Top panel: Obser ved profiles (solid lines and filled dots) of the [Om]A5007 line fitted with two Gaussian components (dashed 
lines) as explained in Section 12,31 Bottom panel: Same as top panel but for galaxies whose profiles are fitted with one 
Gaussian. The plotted errors for the 6dF data (NGC 2369, MGC -03-34-064, NGC 5135, and NGC 7130) are those 
computed from the rms of the continuum adjacent to the line. 

length of this line. While the derived spectr a l reso lution 
is very similar to that from iDasvra et al. I (|2011f) , the 
standard deviation o f our measurements is much lower. 
IDasvra et al. I (|2011|) obtained <7i nst as the average of 
the values of the high excitation fine structure lines of 
a sample of AGN that were deemed not to be resolved 
by Spitzer/IRS with their method. By doing so, they 
probably also included measurements of barely resolved 
lines. This would therefore explain the higher dispersion 
of their instrumental value. 

For the LIRGs we need a criterion to assess whether 
the lines are resolved or not since a significant num- 
ber of sources have measured widths slightly larger than 
the instrumental resolution. Figure Q] shows the distri- 
bution of the measured FWHM of the [Ne in] 15.56 /jm 
line for the sample of LIRGs compared with those of the 
calibration sources. To avoid analyzing barely resolved 
lines we adopt the following criterion for the velocity dis- 
persion to determine whether a line is clearly resolved: 



(cobs - cr inst ) > 3 x v^m + e inst > where tTinst and e inst are 
the instrumental resolution and its standard deviation, 
and dobs and e m are the measured values for the LIRGs. 
The total error of the observations e b s includes both the 
error in the measurement e m and that associated with 
the instrumental value, that is, e obs = (e 2 n + e^t) 1 / 2 - 

Using the criterion described above we found that 
19 LIRGs in our sample show resolved [Ne in] 15.56 /mi 
lines. Most are LIRGs classified as Seyfert (Figure [I]) 
and/or are [Nev] emitters. Those nuclei classified as Hn, 
on the other hand, tend to show FWHM close to the in- 
strumental resolution. Composite nuclei have FWHM 
between those of unresolved lines and the clearly re- 
solved [Ne in] 15.56 /im lines. For those LIRGs deemed to 
have unresolved [Ne in] 15.56 /im lines, the average mea- 
sured FWMH of the [Nelll]15.56 /im line is 0.02776 ± 
0.0017 /im. That is, the observed values for these lines 
are on average 2<r above the instrumental resolution. 

Table [3] lists the intrinsic values of the velocity 
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TABLE 3 

Velocity dispersions and BH masses 



Galaxy 


Class 


ffflKelH]) 


<r([oni]) n 

km s~ 


7([0III]) b 




Ref 


log M B h 






km s _1 


km s 1 


km s _1 




M 


Seyfert nuclei and [NeV] emitters 




oy l 


1 "V\ 4- 1 9 


91 8 4- 99 








O.UU 


NGC 3690 


Sy2 


167 ± 10 






144 ±11 


1 


7.52 


MCG -03-34-064 


Syl 


368 ± 18 


249 ± 10 


656 ± 26 


155 


2 


7.65 


NGC 5135 


Sy2 


138 ±9 


127 ±9 


305 ± 21 


124 ±6 


3 


7.24 


IC 4518W 


Sy2 


156 ± 12 










7.48 


NGC 5990 


Sy2 


264 ± 11 










< 8.45 


NGC 6156 


[Nev] 


160 ± 13 










7.52 


MCG +04-48-002 


[Nev] 


158 ±7 










7.50 


NGC 7130 


Sy2: 


259 ± 16 


130 ±9 


463 ± 32 


147 ±5 


3 


7.55 


NGC 7469 


Syl 


182 ± 13 


145 ±7 


336 ± 17 


152 ± 16 


4 


7.61 


NGC 7591 


Sy2: 


169 ± 17 










7.62 


NGC 7679 


Syl/Sy2 


102 ±7 


176* 




96 


2 


6.77 


Composite and Hll nuclei 


NGC 1614 


Composite 


157 ± 10 






164 ±8 


5 


7.75 


NGC 2369 


Composite 


133 ±8 


146 ± 15 










IC 694 


LINER 


< 70 






141 ± 17 


5 


7.48 


NGC 3256 


Hn 


< 112 






100 ±6 


3 


6.84 


MCG -02-33-098W 


Composite 


103 ±9 












ESO 173-G015 




102 ±7 












IC 4280 


Hll 


128 ±8 












IRAS 17138-1017N 


Composite 


< 80 






72 ±5 


6 


6.24 


IRAS 17578-0400 




119 ±8 












NGC 6701 


Composite 


< 159 






150 ± 20 


7 


7.59 


NGC 7771 


Composite 


178 ± 17 













Notes. — The references for the spectral class of the nuclei are listed in Alonso-Hcrrcr o et al. 1 (|2012D except for 
CGCG 468-002-NED01, which is from the newly analyzed optical spectrum (see Section f2.3[) . 

"The [O hi ] vel ocity dispersion fo r N GC 7679 is also from reference 2. 

References for <x«: 1. IHo et al. I (|2009|L 2. IGu et al. I (120061). 3. IGarcia-Rissmann et al. I ([200511. 4 IQnken et al~l pool ). 5. 
IHinz fc Riekel (|2006I 1. 6. IShier fc Fisher I (|1998I 1. 7. IMarquez et al. 1 (11996,. 



dispersion for the LIRGs with resolved [Ne in] 15. 56 /im 
lines as well as their spectroscopic classes from 
optical spectroscopy and /or mid-IR indicators 
(|Alonso-Herrero et al. I I2012D . The velocity disper- 
sions (corrected for instrumental resolution) are between 
102 and 368kms~ . For galaxies in our sample with 
literature values of the stellar velocity dispersions whose 
[Ne in] 15.56 /zm lines are deemed to be unresolved, we 
provide in this table upper limits to the gas velocity 
dispersion. In Figure [2] we show the normalized profiles 
of the twelve Sy/[Nev] emitter LIRGs with spectrally 
resolved [Ne in] 15.56 /im lines compared with those of 
an unresolved profile represented as a Gaussian function. 

2.3. Optical [Oiu]\5007 measurements 

We obtained archival optical spectra of six LIRGs 
classified as AGN or Composite fr om the six-degree 
Field (6dF) Galaxy Survey (6dfGS, IJones et al. Il200i 
l2009f) and from observ a tions with the FAST sp ectrograph 
(|Tokarz fc Rolil [l997t iFabricant etalH [l998h from the 
Astrophysics (CfA) telescope data center. 

NGC 2369, MCG -03-34-06, NGC 5135, and 
NGC 7130 were observed with the 6dF multi-object fibre 
(angular diameter of 6'.' 7) spectrograph on the United 
Kingdom Schmidt Telescope (UKST). At the typical 
distance of these LIRGs the spectra cover the central 
^2 kpc. The instrument spectral resolution around the 



[Oni] emission lines is a lns t ~125kms 1 . Long-slit spec- 
troscopy was available for CGCG 468-002-NED01 and 
NGC 7469 taken with FAST on the Mt. Hopkins Till- 
inghast 60-inch telescope as part of different observing 
programs. The slit width was 3" and the spectral resolu- 
tion was Cinst ~75kms _1 . None of the spectra are flux 
calibrated, but they can be used to measure the velocity 
dispersion of the emission lines as well as their ratios. 

We first measured the optical line ratios of CGCG 468- 
002-NED01 since we did not have an optical clas- 
sification for this ga laxy. We classify thi s galaxy 
as a Seyfert 2. In iPereira-Santaella et al. I (|2010bl ) 
and lAlonso-Herrero et al. I (|2012D we detected the mid- 
IR [Nev] lines in this galaxy and measured a high 
[0iv]25.89 jum/ [Nell] 12.81 pirn line ratio. We used both 
as evidence of the presence of an AGN in this nucleus. 
The optical classification confirms this. 

For all the galaxies we modeled the [O in] emission 
lines at 5007A and 4959A by fitting a Gaussian profile 
to each line. The relative position of the Gaussians was 
fixed according to the rest-frame wavelengths of the two 
lines. Likewise, the relative intensity was also fixed to 
the value determined by the atomic parameters. We 
found that the fit with one component was only satis- 
factory for CGCG 468-002-NED01 and NGC 2369. For 
the rest we had to add a broad component to the [O in] 
lines. The wavelengths and intensities of the narrow (also 



G 
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called core) component and the broad component were 
considered independent. To estimate the errors of the 
velocity dispersions we used the rms of the continuum 
adjacent to the line for the 6dF spectra and the errors of 
the spectra for the FAST data. The typical errors of the 
velocity dispersions of the [Om] lines are less than 10%, 
and are listed for each galaxy in Table [3J 

The top panel of Figure [3] shows the galaxies for 
which the [O iii]A5007 line fits required two components, 
whereas the bottom panel shows the profiles of NGC 2369 
and CGCG 468-002-NED01. Table [3] lists the velocity 
dispersion of the core of the line and the broad com- 
ponent, if present, corrected for instrumental resolution. 
The velocity dispersions of the core components mea- 
sured in the four galaxies with double components are 
in good agreement with literature er* values, except for 
MCG -03-34-064. 

For the four LIRGs with double components the broad 
component is blueshifted typically by 50 — 300 km s -1 . 
The presence of such asymmetric or broad blue wings 
in the [O ill] A5007 lines is common in AGN (see e.g., 
IGreene fc Hoi [2005] ). The total [O in] emission is dom- 
inated by the broad component in MCG-03-34-064 and 
in NGC 7130, whereas in NGC 5135 and NGC 7469 the 
core component domin ates the emission. In NGC 5135 
iBedregal et al. I (|2009f ) based on both the relatively 
broad [Si Vi] 1.96 /mi coronal line and its spatial extent 
suggested the existence of different kinematical compo- 
nents and the possible presence of A GN induced out- 
flows. Gonzalez Delgad o et al. I ([1998D found blueshifted 
components in the Lya and [Oni]A5007 lines of both 
NGC 7130 and NGC 5135. The velocity dispersion of the 
broad component of the [O in] A5007 line of NGC 7130 is 
similar to the Ha bro ad component (er ~ 400 km s -1 , 
IBellocchi et al~l [2011 . In NGC 7469 IWilson et al.1 
(|1986f ) also detected two components in the [Oiii]A5007 
line. 

3. RESULTS 
3.1. Black Hole Masses 

In Tabled we list, in addition to the gas velocity disper- 
sions from the [Nem]15.56 /mi and optical [Oiii]A5007 
lines, the stellar velocity dispersion cr* values from the lit- 
erature, when available. We also included in this table cr* 
values from the literature for other LIRGs in our sample 
(IC 694, NGC 3256, IRAS 17138-1017, and NGC 6701) 
whose [Ne in] 15.56 /im lines appear unresolved at the SH 
spectral resolution. We note, however, that in the case 
of on-going major mergers the stellar velocity dispersion 
v s. BH mass rel a tion might n ot be applicable. 

IDasvra et al~l (|2008l 120 111 ) showed that mid-IR re- 
solved fine-structure lines can be used to study the NLR 
of AGN. In particular they demonstrated that the mass 
of the BH correlates with both the velocity dispersion 
and the luminosity of the NLR region for a sample of 
optically selected AGN. The relation is: 

log (Mbh/M ) = -1.82 + 4.24 log( ( r NcIII /kms- 1 ) (1) 

IDasvra et al. I (|2011[ ) fitted a similar relation for the 
optical [Oni]A5007 line: 



log(M BH /M ) = -1.78 + 4.24 log^om/kms- 1 ) (2) 

For both relations IDasvra et al. I (|2011l ) fixed the 
slope to the value derived from the correlation between 
the stellar velocity d ispersion and the BH mass by 
IGultekin et al~l ((20091) : 

log (Mbh/M ) = 8.12 + 4.24 log(cr >t /200kms" 1 ) (3) 

The typical uncertainties on the derived Mbh can be as 
high as 0.8 de x if using the gas a based on the rms of 
the relations ( | Dasvra et al. 1 1201 lT ) and 0.4 dex if using 
a* (|Gultekin et al. II2009D . ~ 

Before we compute the BH masses, we compare 
the gas velocity dispersions for the LIRGs with both 
[Ne in] 15.56 /im and [O in] A5007 measurements. As we 
saw in Section 12.31 in some cases the optical line is fit- 
ted with two components. For AGN there is generally a 
good agreement between cr* and the gas v elocity disper- 
sion from the narrow comp onent of [O in] (|Onken et al. I 
l200i IGreene fc Ho I [20051) . Therefore, the core of the 
[Om]A5007 line probes the gravitational movements in 
the NLR and thus can be used to estimate Mbh- On the 
other hand, the broad components tend to be blueshifted, 
may be more influenced by the AGN, and are likely pro- 
duced by outflows. As can be seen from Tabled the ve- 
locity dispersions from the [Ne in] 15. 56 /im line appear 
to be intermediate between those of the two optical com- 
ponents. Spectrally resolved blueshifted [Ne in] 15. 56 /mi 
lines (as well as th e [Ne v] lines) have al so been reported 
for local ULIRGs ([Spoon fe Holt 1 120091) . and have been 
interpreted as the result of outflows. Thus, it is possible 
that in some cases we might overestimate the BH masses 
when using the [Ne in] 15. 56 /mi line if there are several 
components. The most suspect galaxy in our sample 
would be NGC 5990, for which (TNem predicts the most 
massive BH by far in our sample. We note, however, that 
the calibration of Dasyra et al. (2011) should take these 
effects into consideration. 

In Table [3] we list Mbh as calculated from, in this 
order of preference, the stellar velocity dispersion, the 
velocity dispersion of the core of the [O in] line, and 
the [Ne in] 15.56 /xm line velocity dispersion. For those 
LIRGs with no definitive evidence of AGN activity 
(i.e., classified as composites, Hn or unknown), we 
only list Mbh in the second part of the table if we 
have a value of <7*. This is because we cannot be 
sure whether the [Neni]15.56 /mi emission comes mostly 
from the NLR or whether it is prod uced by star forma- 
tion (|Pereira-Santaella et aT~ll2010bl ). 

For those LIRGs hosting an actively accreting BH (that 
is, with clear signs of AGN activity, first part of Table [3]) 
we find typical BH masses of 3 x 1O 7 M , ranging from 
6 x 1O 6 M to 3 x 1O 8 M . The estimated BH masses 
for those LIRGs classified as composite or Hll are also 
in this range. The masses of the LIRG BHs appear to be 
similar to those of the currently growin g BH in the lo- 
cal universe hosted in late- typ e galaxies (|Heckman et al. I 
2004; ISchawinski etlffll2l)l"(ih . 

The typical Eddington ratios (Lboi(AGN)/LEdd) f° r 
the LIRGs with an estimate of Mbh are 2 x 10~ 2 , and 
range between 5 x 10~ 4 and 0.1. This means that the 
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BH in local LIRGs are accreting at a lower efficiency 
than those in local ULIRG (typical values of 0.08 — 0.4 
depen ding on the method used for estimating the BH 
mass, IDasvra et al~l 120061 : IVeilleux et al. 1120091 ) . 

The majority of LIRGs with spectroscopically resolved 
[Ne in] 15.56 /im lines are classified as Seyfert galaxies. 
However, in local LIRGs there is a high fraction of com- 
posite objects (AGN/starburst). Moreover, this fraction 
remains approximately constant as a function of IR lumi- 
nosity from galaxie s with £ir ~ 10 10 £© up to ULIRGs 
([Yuan et al. 1 12010). As can be seen from Figure[TJ com- 
posite nuclei in local LIRGs have values of the FWHM of 
[Ne in] 15.56 /urn on average lower than those of Seyfert 
nuclei but higher than those of nuclei classified as H n 
and the instrumental resolution. Assuming that these 
composite nuclei do indeed host an AGN and that the 
[Ne in] 15.56 /im line is probing their NLR, then their 
BH masses ought to be less than ~2-3x 10 7 M©. The 
estimated typical AGN bolometric luminosities from the 
mid-IR spectral decomposition are 10 43 ergs -1 or less, so 
we cannot constrain their Eddington ratios. 

The BH of local LIRGs are only marginally less massive 
than those of local ULIRGs. The B H masses of the latter 
are in the range ~ 10 7 - 5 x 10 s (iTacconi et al. 1120021: 



tends typicall y 1.5 kpc. We use the following recipes 
put forward bv lDiamond-Stanic fc Riekel (|2012l ) to com- 
pute the nuclear SFRs based on the [Nen]12.81 /im line 
11.3 /im PAH feature luminosities: 



IDasvra et al. I2006t IVeilleux et al. 120091) . Samples of lo- 
cal ULIRGs are mostly domi nated by gas-rich interactin g 
galaxies and major mergers USanders fc Mirabel I [[996) , 
with coalesced ULIRGs havin g slightly larger BH m asses 
than pre-coalescence ULIRGs (|Dasvra et al. 1 12006). Our 
sample is both flux and volume-limited and thus is com- 
posed mostly of LIRGs with £ir ~ 1 — 2 x 10 11 £© (in- 
dividual galaxies), which also tend to be spiral galaxies, 
minor mergers, and galaxies in groups. That is, morpho- 
logically local LIRGs are not dominated in numbers by 
major mergers (see e.g. [S anders fe Ishida ll2004tlKavirai I 
20091 IPereira-Santaella 112011 " and their BH masses are 
simil ar to those of pre-coalescent ULIRGs (Da svra et al. I 
[20061) . 

3.2. Star Formation Rates 

We computed the SFRs for our sample of LIRGs 
using a number of IR-based indicato rs, including 
the fi ne-structure [Ncn] 12.81 /im line (|Roche et al. I 
[19911: iHo fc Keto I l2007n. the 11.3 /im PAH featur e 
(|Brandl et al. I 120061: iDiamond-Stanic fe Riekel I2012D 
and the total 8 — 1000 /xm IR luminosity ([Kennicuttl 
1998). Although it is not clear whether the presence 
of an AGN could destroy the PAH carriers or not, at 
the typical luminosities of the AGN hosted in local 
LIRGs the 11.3 PAHs do n ot seem to be affected 
(|Diamond-Stanic fe Rie"ke~1l2010D . 

For the integrated SFRs we used £ir of the indi vidual 
galaxies as given in lAlonso-Herrero et al. I (|2012t their 
table 1) after subtr a cting the AGN contribution. We 
used the iKennicuttl ([1998D relation in terms of the IR 
luminosity converted to a Kroupa IMF: 

SFR(M© yr' 1 ) = 1.14 x KT 10 (£ m , £©). (4) 

We obtained total SFRs for the individual galaxies in our 
sample of between 1 and 60M©yr . 

For the nuclear SFR we used the Spitzer/TRS spec- 
troscopy. At the median distance of 65 Mpc for our sam- 
ple of LIRGs, the IRS SH slit width (4.7arcsec) sub- 



SFR(M© yr" 1 ) = 9.6 x 10~ 9 £(11.3 /mi, £ c . 



(5) 



SFR(M© yr" 1 ) = 8.9 x 10" 8 £([NeII], £©). (6) 



Diamon d-Stanic fc Riekel (|2012[ ) calibrated Equa- 
tio ns El and [H fo r gala xies with £ir < 10 11 £© using 
the iRieke et al. I (|2009D templates and a Kroupa IMF. 
We chose to use these calibrations rather than those for 
higher IR luminosities because the median value of the 
IR luminosity of the individual galaxies in our sample is 
1.3 x 10 11 £©. In addition, the width of the IRS SH slit 
only probes the nuclear regions of our sample of LIRGs 
and thus lower IR luminosities. 

We also note that Equation [5] assumed that the flux 
of th e 11.3 iim PAH fea ture was measured using PAH- 
FIT ([Smith, et al. I 120071) . The spectral coverage (~ 
10 — 18 /im) and spectral resolution of the SH data used 
here are not adequate to use pahfit. We measured 
the PAH fluxes usin g a local conti n uum (Section [2] and 
Table g]). Following [Smith et al. I ([2007D we applied a 
multiplicative factor of two to the PAH fluxes measured 
with a local continuum to make a proper comparison 
with pahfit fluxes. We also corrected the 11.3 /im PAH 
fluxes for extinction using the nuclear strength of the 
10 /im silicate feature (Lrom lAlonso-Herrero et al. 1120121) 
and the extinction law given in ISmith et al. I (120071) . 
For most of our local LIRGs this correction is small, 
as the s ilicate absorptions are mode rate (see e.g. fig- 
ure 6 in lAlonso-Herrero et al. | |2012) as compared with 
the more deeply em bedded population of local ULIRGs 
(jSpoon et al. Il2007h . 

While the [Neil] 12.81 /im emission in low luminos- 
ity and moderate luminosity AGN is mostly produced 
by star formation, in high luminosity AGN this emis- 
sion line ca n also have an importan t contri bution from 
the AGN ([Pereira-Santaella et al~l l2010bO . To cor- 
rect the [Nell] 12.81 /im line flux for any possible AGN 
contribution, we follo wed the method put forward by 
IVeilleux et al. I (|2009t ). This method reproduces the ob- 
served values of the [Oiv]25.89 jnm/[Neiill2 .81 /xm line 
ratio (given in lAlonso-Herrero et al. I [20 12f ) as a frac- 
tional combination of the typical line ratio of an AGN 
(~ 4) and that of a starburst galaxy (~ 0.01), where 
the AGN contribution is the free parameter. Using this 
method we find that the AGN contribution to the nuclear 
[Nen]12.81 /im emission of local LIRGs is small. These 
corrections are only necessary for 10 LIRGs and the AGN 
contribution to the nuclear [Neil] 12.81 /im fluxes ranges 
from 3% to 53% (average ~ 18%). 

Figure [4] (left panel) compares the nuclear SFRs as 
derived from the [Neil] 12.81 /im line (after removing 
the AGN contribution) and from the 11.3 /im PAH fea- 
ture. We have color coded symbols by the AGN bolo- 
metric contribution to IR lu minosity of the galaxies (see 
lAlonso-Herrero et al. II2012D . As can be seen from this 
figure, the majority of LIRGs with small or no (<25%) 
AGN contribution show a good agreement between the 
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Nuclear SFR([NeIl]12.B/xm) M yr 
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Fig. 4. — Left panel: comparison between the nuclear SFRs from the [Ncn]12.81 (im line (corrected for AGN emission) and the 11.3/im 
PAH feature for the sample of local LIRGs. The colors denote the AGN bolometric contribution to Lir of the individual galaxies: < 1% and 
no AGN contribution (blue), 1 — 7% (yellow), 7 — 25% (orange), and > 25% (red). The highly discrepant point with no AGN contribution 
to the left of the straight line is IC 860. The straight solid line is a 1:1 line and is the average value of the ratio between the nuclear SFRs 
computed with the two estimators (see the text). Right panel: comparison between the nuclear SFR from the 11.3 ^m PAH featu re and 
the integrated SFR from the IR luminosity after subtracting the AGN contribution as estimated by Alonso-Hcrrcro ct al. (2012). The 
straight solid line shows an average ratio between the nuclear and the total SFRs of 0.5 (see text), typical of local LIRGs, whereas the 
dashed line represents a 1:1 relation. 



derived nuclear SFRs using the two indicators. Indeed, 
we find that the ratio between the nuclear SFR esti- 
mated from the 11.3 /im PAH feature and that from the 
[Nen]12.81 /zm line is on average 1.0 ± 0.3. 

In a few cases with a significant AGN bolomet- 
ric contribution the nuclear SFR derived from the 
[Nen]12.81 /im line is higher than that from the 11.3 /jm 
PAH feature, even after correcting for the AGN contribu- 
tion. Given the good agreement between the nuclear SFR 
computed with the two methods for LIRGs with small or 
no AGN contribution and to avoid AGN contamination 
issues from now on for the nuclear SFR we use those 
derived from the 11.3 /im PAH feature. Thus, the indi- 
vidual nuclei of the sample of local LIRGs show values 
of the nuclear SFR of between ~ 0.8 and ~ 20 M Q yr -1 . 

In Figure @] (right) we compare the total SFR against 
the nuclear (typically on scales of 1.5 kpc) SFR for the 
individual galaxies of the LIRG sample. It is clear 
that a significant fraction of the individual galaxies in 
our sample have more than half of their total SFRs 
taking place outside the nuclear (1 — 2 kpc) regions. 
For the individual galaxies of the sample the average 
value of the nuclear SFR over the total SFR ratio is 
0.5 ± .3. This agrees with findings using other indi- 
cators (lHattori et al. 112004 lAlonso-Herrero et al. Il2006t 



iDiaz-Santos et al. ll2010HRodriguez-Zaurm et al. 11201*1*1 ) 

We also note, however, that in approximately 20% of our 
volume-limited sample of LIRGs the nuclear emission ac- 
counts for most of the total SFR measured in the indi- 
vidual galaxies, as is the case for a large fraction of local 
ULIRGs. 

Incidentally, the same comparison for the nuclear SFRs 
(Figure IH left) but usi ng the SFR equations calibrate d 
for Lir > 10 11 Lq fsee lDiamond-Stanic fc Rieke II2012D . 
results in nuclear SFRs from the 11.3 /im PAH feature 
which are on average a factor of 1.3 higher than those 
from the [Nen]12.81 /im line. Moreover, the ratio be- 
tween the nuclear SFR from the 11.3 /im PAH feature 



and the total SFR from Lir would also be on average 
1.3. This gives us confidence that we used the appropri- 
ate SFR calibrations for the nuclear rates. 

3.3. Relation between BH accretion rates and SFR in 

LIRGs 

The presence of an AGN is an unambiguous signpost 
of a period of BH growth. The AGN luminosity can then 
be expressed in terms of the BHA R ?71rh and the mass- 
energ y conversion efficiency e (see Alex ander fc Hickox I 
2012, and references therein): 

m BH (M yr" 1 ) = 0.15(0.1/e)(i bol /10 45 ergs" 1 ) (7) 

We used this relation wit h the typical value o f e = 
0.1 in the local universe (M arconi et al. I 120041 ) and 
Lboi(AGN) from the mid-IR spectral decomposition (see 
lAlonso-Herrero et al. 112012ft . We obtained BHAR be- 
tween 0.0007 and 0.08 M Q yr- 1 for local LIRGs. The 
uncertainties of the estimates of the BHAR are domi- 
nated by the uncertainties of the AGN bolometric lu- 
minosities which are typic ally less than 0.4dex (sec 
lAlonso-Herrero et al. 1120121 for details). 

Having calculated the BHAR we can now compare 
them with the nuclear and total SFR for the sample of 
LIRGs. The distribution of SFR/tobh can provide clues 
as to whether AGN and on-going star formation activity 
are contempor aneous in l ocal L IRGs. Recent numerical 
simulations by iHopkins I (|2012f) predict a time offset be- 
tween the peaks of these activities which is thought to 
depend on the physical scale within the galaxy and the 
dynamical time of the galaxy, among other parameters. 

Figure [3] shows the distribution of SFR/mRH for the 
sample of local LIRGs for the nuclear (~ 1.5 kpc) and the 
integrated SFRs, only for those LIRGs with an estimate 
of iboi(AGN). The nuclear and integrated SFR/m B R 
are in the range ~ 12 to ~ 2 x 10 4 , with median values of 
the nuclear and integrated log(SFR/mBH) of 3.1 and 3.4, 
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Nuclear (~1.5kpc) log(SFR/m m ) 

Fig. 5. — Distributions of the SFR/riiBH ratios for nuclear (typi- 
cal scales of 1.5 kpc) SFRs (lower panel) and integrated SFRs (up- 
per panel) for the sample of local LIRGs with an estimate of the 
AGN luminosity. The nuclear SFR are from the 11.3 ^tm PAH fea- 
ture and the integrated SFR from Ljr after subtracting the AGN 
component (see the text for details). 
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IR-bright in RSA 
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Fig. 6. — Distributions of the ratio between the nuclear SFR and 
the BHAR (bottom panel) and the integrated SFR and the BHAR 
(upper panel) for the RSA sample of AGN (blue). In orange we 
mark those RSA AGN that are also IR-bright galaxies (see text). 
We only included RSA galaxies with ttibh > 10 — 4 A^q yr _1 . 



respectively. Most of the nuclei classified as Seyfert from 
optical spectroscopy show nuclear log(SFR/m,BH) < 3.0. 
The nuclei classified as composite tend to show larger 
values of the nuclear SFR to BHAR ratio, as expected if 
they were powered by both SB and AGN activity. 

For the entire sample of local LIRGs we obtain 
5^SFR/^mBH > 2900. This is a lower limit because 
for half of the sample we did not detect an AGN com- 
ponent and therefore we are only using an upper limit 
for their BHAR. If we repeat this only for those LIRGs 
with an estimate of their AGN bolometric luminosities 
we obtain a time- averaged value of SFR/mBH — 2500, 
if we assume that all LIRGs go through an AGN phase, 
that is, a ~ 50% duty cycle. The two values are con- 
sistent with each other taking into account the typi- 
cal uncertainties in calculat ing L\, \{AGN) (~ 0.4 dex, 
lAlonso-Herrero et al. I [20121. This ratio is a few times 
higher than that measured for t he local population o f 
bulge-dominated galaxies (~ 10 3 , iHeckman et al. Il2004t ) 
and the local no rmalization of the BH mass versus bulge 
mass relation ([Marconi fc Hunt 1 l2003t lHaring fc Rix I 
12004 . Even though the SFRs of LIRGs are much higher 
than those of local bulge galaxies, the time-averaged 
SFR/mBH value of local LIRGs is similar to that of local 
bulge gala xies during the fir st 0.3 Gyr of their on-going 
starburst ([Wild et al. 1 120101 their figure 9). This sug- 
gests the existence of similar time delays between the 
peak of star formation and BH growth in local LIRGs as 
well. 

4. DISCUSSION 

4.1. Comparison between the AGNs in LIRGs and the 
optically identified RSA Seyferts 



In this section we compare the properties of the AGNs 
identified in local LIRGs with the optically selected 
Seyfert galaxies in the revis ed-Sha pley-Ames catalog, 
(RSA , iSandage fc Tammannl 119871 : iMaiolino fc Rieke I 
1995). This optical sample was spectroscopically selected 
from the original galaxy-magnitude limited sample of 
RSA galaxies and is believed not to be biased against 
low-luminosity Seyfert galaxies. Also, as for the sample 
of local LIRGs studied in this work, the Seyfert activity 
in the RSA sample does not appear to be driven primar- 
ily by mergers for the majority of the galaxies. Finally, 
the AGN bolometric luminosities of the RSA Seyferts are 
similar to those of the local LIRGs. 

The median values of SFR/mBH for our LIRGs are 
on average higher than those of local Seyfert galax- 
ies, especially the nuclear va lues (see Figure 5 and 
IDiamond-Stanic fc Rieke 1 1201 21. Of course, this is not 
completely unexpected as local LIRGs, including those 
hosting an AGN, are essentially selected by their star for- 
mation activity, whereas the RSA Seyferts are selected by 
their host galaxy magnitudes. The latter is more likely 
probing the stellar mass, rather than the star formation 
activity. 

We can make a more meaningful comparison if we iden- 
tify the IR-bright Seyferts in the RSA sample. We set 
the limit to log(LiR/LQ) = 10.8 to match approximately 
the lower limit covered by the IR luminosities of the indi- 
vi dual galaxies of our sample o f local LIRGs (see table 1 
in lAlonso-Herrero et al. Il2012| ). We also included from 
our sample the newly identified Seyfert nuclei with an 
estimate of Lboi(AGN). The values of SFR/ msH for the 
RSA Seyferts are taken from IDiamond-Stanic &: Riekel 
(2012) except for those already included in our sam- 
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pie of LIRGs. iDiamond-Stanic fc Riekel ([2012) com- 
puted the nuclear SFRs from the 11.3 PAH fea- 
ture and the MIPS 24 /im photometry, and L bo i(AGN) 
from the [0iv]25.89 /im line. Finally, for the few RSA 
Seyferts deemed to have their [Oiv]25.89 /im luminosi- 
ties strongly contaminated by star formation, we esti- 
mated Lboi(AGN) (and thus BHAR) from their hard 
X-ray luminosities and using a bolometric correction. 
In this comparison we only included RSA Seyferts with 
^bh > 10 _4 Af©yr _1 , to match the approximate AGN 
detection threshold in the LIRG sample. With this 
limit the IR-bright galaxy fraction in the RSA sample 
of Seyfert galaxies is then ~ 22% ± 5%. 

Figure [6] shows the distributions of the nuclear and 
total SFR/meH for the all the RSA Seyferts indicating 
which ones are also IR-bright. Clearly, those RSA Seyfert 
classified as IR-bright have ratios similar to, although 
slightly smaller than, those of the complete sample of lo- 
cal LIRGs. Given the fact that the AGNs in the RSA 
sample are not very different from those in the LIRGs, 
the results above would suggest that the bright AGN 
phase comes after and is somewhat distinct from the 
LIRG star forming phase. This interpretation would be 
consistent with observational works showing a delay be- 
tween the onset of the star formation activity and the lat- 
ter feeding of the AGN with the conseq uent BH growth 
(|Davies et al. 1120071 iWild et al. 1120101) as we ll as with 
predi ctions from numerical simulations (e.g., IHopkinsI 
l20ll . 

4.2. SFRs versus AGN Detection in local LIRGs 

It is well known that the AGN fraction increases 
with increasing IR luminosity from LIRGs to ULIRGs, 
and at the highest IR luminosities the AGN might 
domi nate bolometric a lly t he lum i nosity of th e sys- 
tem (IVeilleux et al I [1991 [200l lYuan et al. I l2010l: 
iNardini et al. 1 120101: lAlonso-Herrero et al. 1 120121 ) . We 
now examine further the possibility that the IR-bright 
AGN phase comes after the distinct IR-bright star- 
forming phase. If that were the case, we would expect 
the AGN fraction in LIRGs to be higher at the lowest 
SFR. The AGN fraction here is computed only for those 
LIRGs with a secure AGN detection, that is, a Seyfert 
classification and/or the presence of mid-IR [Nev] lines. 
As in the previous section, we computed the integrated 
SFR of the individual galaxies from Ltr after subtracting 
the AGN contribution. 

Figure [7] shows the AGN fraction as a function of 
the total SFR of the galaxy. There is a tendency for 
LIRGs in lowest SFR bin to have a higher AGN in- 
cidence (36%^2%) than those in the highest SFR bin 
(17%^°%). Alternatively, this could be explained if it 
is harder to detect AGN in LIRGs with very high SFRs, 
especially if it takes place mostly in the nuclear regions. 
Conversely, it may be easier to identify an AGN for low 
nuclear SFRs and extended star formation. The first ex- 
planation is not likely however, as most of the LIRGs in 
our sample in the high SFR bin also tend to have rather 
extended SF. In other words, from Figure 0] (right panel) 
there is no tendency for the fraction of nuclear SFR to 
increase with the total SFR of the galaxy. In fact, in 
local LIRGs the mid-IR emission (both from AGN and 
star formation activity) only starts to be highly concen- 




20 30 40 50 
Total SFR(L IR ) M e yr" 1 

Fig. 7. — Distributions of the integrated SFRs of the local LIRGs 
(grey histograms) and of those with a Seyfert classification and/or 
[Nev] emitters (solid line histograms). The data points represent 
the AGN fraction (scale on the right hand vertical axis) in SFR 
intervals with la error bars for the AGN fraction. 



trated in the central re gions for IR luminosities above 
10 n - 8 L Q , as shown by IPfaz-Santos et al. I ([2010). We 
conclude that the higher AGN incidence at low SFRs in 
local LIRGs provides further evidence for a time delay 
between the peaks of the star formation and BH growth. 

4.3. Key differences between merger and non-merger 
local LIRGs 

Before we explore the role of the LIRG phase in the 
context of star formation activity and BH growth of mas- 
sive late- type galaxies (next section), it is important to 
point out the differences in the star formation histories 
of merger and non-merger LIRGs in the local universe. 

The triggering of the activity in local ULIRGs and 
likely also in the most luminous local LIRGs is driven 
by major mergers. As a consequence, ULIRGs tend to 
host more luminous AGN (quasar or nearly quasar-like 
luminosities), slightly more massive BH and have higher 
AGN bolome tric contributions that local LIRGs (see Sec- 
tion IO and IVeilleux et al. I [200l I Nardini et al~l I20T0I: 
Alonso-Hcrrcr c~et al. II2012D . Local merger LIRGs do 
show evidence for a recent and intense period of star for- 
mation that c onsumed the gas faster, that is, a bursty 
SFR (see e.g., lAlonso-Herrero et al. M2000L [20 01) as in 
local ULIRGs (|Rodriguez-Zaurin et al. I I2010D . In this 
respect, most of the local merger-LIRGs could be consid- 
ered as sub- ULIRGs and they might go through a future 
ULIR G phase or may alrea dy have experienced such a 
phase ([Murphy et al. 1120011 ) . 

Major mergers, however, do not dominate in num- 
bers the population of local LIRGs, especially at Ltr < 
3 x 10 11 Lp), as inferred from morphological studies 
(see e.g.. iSanders fc Ishida 1 120041 : lAlonso-Herrero et al. I 
120061 : iKavirai 1 120091 ) . Therefore, in most local LIRGs 
the current episodes of star formation activity and BH 
growth are likely the result of a less violent process, 
e.g., via minor mergers, fly-by companion s, and/or sec- 
ular evolution. iPoggianti fe Wu I ([2000D showed that 
most isolated galaxies in their sample of IR-bright galax- 
ies showed on average more moderate optical Balmcr 
absorption features than the strongly interacting sys- 
tems in their sample. This suggests that in non-merger 
LIRGs there is not a strong contribution from a dominant 
post-starburst stellar population that would be associ- 
ated with a recent bursty star formation history. This 
agrees with the finding that most of the mass in non- 
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Fig. 8. — Cartoon (not drawn to scale) depicting the proposed 
scenario the current episode of star formation and BH growth of 
a local non-major merger LIRG with Lj^ < 3 X 10 11 Lq. The 
last episode of SF is assumed to be an exponentially decaying SFR 
(red curve) with an e-folding time of r = 2.5 X 10 s yr. The current 
period of BH growth is assumed to have a 5 X 10 8 yr lifetime (blue 
curve), as in local galaxies with A/bh < 10 s Mq ((Marconi et al. I 
120041 1. The BHAR would be on average ~ 3000 times lower than 
the SFR during the IR-bright phase. 

merger LIRGs is from evolved stellar populations formed 
a few Gyr ago with a much smaller contri bution in mass 
from a relatively young stellar pop ulation (|Kavirai Il2009t 
lAlonso-Herrero et al. 1120061 [2010 ) . 

4.4. The role of the LIRG-phase in the growth of BH at 
Lib. < 3 x 10 11 L Q 

The cartoon in Figure [8] illustrates our proposed sce- 
nario for the current episode of BH growth in a non- 
merger local LIRG typically with Lir < 3 x 10 11 Lq 
(i.e., integrated SFR< SOMeyr^ 1 ). This IR lumi- 
nosity marks approximately the border between lo- 
cal populations of LIRGs dominated by major merger 
and populations of disky and relatively isolated LIRGs 
(ISanders fc Ishida ll200l . 

We assume that the galaxy has an old stellar 
population formed a few Gyr ago (jKavirai I 120091 : 
lAlonso-Herrero et al. I [2010h and an existing BH. Since 
the galaxy did not experience a recent major merger 
but rather a less violent process that induced the last 
episode of SF, we can assume an exponentially de- 
caying SFR for it, as SFR(i) = SFR e~*/ r where 
r is the e-folding time. The peak SFR(£ = 0) ~ 
30M Q yr _1 is driven by the Lir limit set above. From 
the p osition of local AGN on t he color-magnitude dia- 
gram, iSchawinski et al.~1 (|2009f) demonstrated that the 
e-folding time cannot be too short (a few Myr or bursty) 
or too long (i.e., const ant SFR) . The high AGN detection 
rate i n local LIRGs dVeilleux et al"~l Il995t lYuan et al. I 
120101: lAlonso-Herrero et al. I I2012D implies that LIRGs 
are able to feed the BH (i.e., show AGN activity) roughly 
during half of the LIRG-phase (defined as having in- 
tegrated SFR> 10A/ Q yr _1 , see Figure S]). Moreover, 
given the moderate BH masses inferred for LIRGs (Sec- 
tion [3Tj, this AGN activity is likely to be sustained for 
~ 5 x 10 s yr, that is, the typical lifetime of the BH growth 
for M BH < 10 8 M Q (see lMarconi et al. Il200l . 

In this simple scenario, the galaxy would be first clas- 



sified optically as an H II LIRG and would show elevated 
integrated and nuclear SFRs. Provided that there is a 
mechanism to bring the gas near the BH gravitational po- 
tential, after a time the existing BH starts accreting ma- 
terial at a sufficiently high rate. Then the galaxy would 
still have a Seyfert-like AGN luminosity. Because there 
is an elevated nuclear SFR, the galaxy would be optically 
classified first as composite (AGN/SB) and later on as a 
Seyfert LIRG as the SFR decreases. In both cases the 
SFR/meH ratios would be on average higher than those 
of optically identified Seyferts. After a few hundred mil- 
lion years the SFR would be below 10 M Q yr -1 and thus 
it would not be classified as a LIRG any longer. However, 
the AGN phase may be still active for another few hun- 
dred million years as Seyfert-like AGN s do not requir e 
large amounts of gas to remain active (|Hopkins II2012T ). 
At this point the galaxy would be identified as an optical 
Seyfert with a lower value of the SFR/BHAR ratio and 
a moderate SFR. Eventually both the SF and the BHAR 
would be low and the galaxy may not be identified as a 
Seyfert. A SFR e-folding time of r ~ 2.5 x 10 8 yr and a 
BH growth lifetime of ~ 5 x 10 8 yr would reproduce the 
SFRs as well as the Seyfert an d composite fractions of 
local LIRGs (jYuan et al. \MHh . 

In summary, in the local universe LIRGs with Lir < 
3 x 10 11 Lq LIRGs and an AGN represent an early phase 
of the (possibly episodic) growth of BHs in massive spi- 
ral galaxies with high SFR, not necessarily associated 
with a major merger event. An Hn-like LIRG phase 
would predate the AGN classified LIRG phase. The 
latter in turn would be followed by an optically iden- 
tified AGN phase where the galaxy still shows an ap- 
preciable SFR but below the LIRG limit. Although 
the triggering mechanisms for the SB and AGN activ- 
ity of LIRGs and ULIRGs might be different, it ap- 
pears that the sequence of events is similar. Namely, 
ULIRGs are believed to involve the interaction between 
two gas-rich galaxies that results fi rst in a high SFR and 
then an optical QSO phase (see ISanders et al. 1 119881 : 
ISanders fc Mirabel 1 119961 : lYuan et al. 1 120101 and refer- 
ences therein). The end product, in contrast, is rather 
different with U LIRGs probably evolving into moderate- 
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5. SUMMARY 

We have studied the current co-evolution of the star 
formation activity and BH growth for a complete sample 
of 51 local LIRGs (see lAlonso-Herrero et al~1l2006l [20121 
for details on the sample) . We used Spitzer /IRS SH (R ~ 
600) data in the ~ 10 — 18 /im spectral range to estimate 
BH masses from resolved [Nein]15.56 /im fine structure 
lines, and to measure nuclear SFRs from the 11.3 /zm 
PAH feature and the [Neil] 12.81 /im line. 

We detected spectrally resolved [Nein] 15.56 /im lines 
in 19 LIRGs with gas velocity dispersions between 102 
and 368 km s -1 . Most of these are LIRGs classified as 
Seyferts and/or are [Nev] emitters. We assumed that 
[Ne in] 15.56 /im probes gas ari sing from the NLR in the 
AGN LIRGs and then used the lDasvra et al. I (|2011h re- 
lation between Mrh and [Ne in] 15. 56 /im velocity dis- 
persion. We combined gas velocity dispersions from 
[Ne III] 15. 56 jum and [Oiii]A5007A and stellar velocity 
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dispersions from the literature to obtain BH masses be- 
tween M B h = 6 x 10 6 M Q and M B h = 3.5 x 10 8 M (me- 
dian value of 3 x 10 7 M Q ). The derived AGN bolometric 
luminosities of L bol (AGN) = (0.4 -50) xlO 43 ergs" 1 (me- 
dian of ~ 1.4 x 10 43 ergs _1 , see lAlonso-Herrero et al. I 
IMl imply BHARof m B H ^ 0.0007-0.08 A^yr" 1 and 
relatively low Eddington ratios (typically 2 x 10~ 2 ). 

Although AGNs accompany the star f ormation activ- 
ity in the large majority o f local LIRGs (IVeilleux et al. I 
fl995t lYuan et al. I Mm I Alonso-Herrero et al. I I2012h~ 7 
in most cases the AGN is n ot energetically dominant 
(|Alonso-Herrero et al. II2012D . Using the 11.3 /zm PAH 
feature we obtained nuclear (on typical scales of 1.5 kpc) 
SFRs for the individual galaxies of between ~ 0.8 and 
~ 20 Mq yr _1 and using Lir integrated SFRs between 1 
and 60M Q yr _1 (for a Kroupa IMF). In a large frac- 
tion of local LIRGs this on-going star formation ac- 
tivity is not only taking place in the nuclear regions 
but it is also spread out on kpc scales throughout the 
galaxy. From the IR estimates in this work the median 
value for the sample of local LIRGs is SFR(nuclear ~ 
1.5kpc)/SFR(total)= 0.5. 

We found that local LIRGs have a time-averaged value 
of SFR/mBH > 2900, which is a few times higher than 
that of the l ocal population of bul ge-dominated galax- 
ies (~ 10 3 , Heck man et al. 1 120041 ) and the local nor- 
malization of the BH mass versus bul ge ma ss relation 
(jMarconi fc Hunt 112001 lHaring fc Rix 112004ft . However, 
this ratio is similar to the value of local bulge galaxies 
averaged during th e first ~ 0.3 Gyr of the last starburst 
()Wild et al. 1120101 ). that is, before the peak of the BH 
growth. 

Local LIRGs have on average higher SFR/ttibh ratios 
(both nuclear and total SFR) than the optically selected 
RSA Seyferts in the local Universe. However, the IR- 
bright RSA Seyferts have similar ratios to those of local 
LIRGs. Since the overall properties of the AGNs in local 
LIRGs and in the RSA sample are very similar, we con- 
cluded that the AGN phase comes after and is somewhat 
distinct from the LIRG star forming phase. The AGN 
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fraction in local LIRGs appear to be a function of the in- 
tegrated SFR, with LIRGs in the low SFR bin having a 
higher AGN (i.e., Seyfert) incidence than those with high 
SFRs. This provides further support to the scenario that 
in local LIRGs the peak of intense IR-bright SF activity 
occurs prior to the AGN phase. 

Finally we put forward a simple model to explore the 
role of the LIRG-phase in the last episode of BH growth 
at Z/ir < 3 x 10 11 Lq. At these luminosities most LIRGs 
are not major-mergers, and have an old stellar popula- 
tion formed a few Gyr ago and a current episode of star 
formation. The latter only a small mass fraction of the 
total stellar mass. We propose that an exponentially de- 
caying SFR with a peak SFR(i = 0) ~ SOMgyr" 1 and 
an e-folding time of r ~ 2.5 x 10 8 yr together with a 
5 x 10 8 yr period of BH growth would explain the ob- 
served SFRs and fractions of Seyfert and composite ob- 
jects in local LIRGs. In this scenario non-merger LIRGs 
hosting an AGN represent an early phase of the possibly 
episodic growth of the BH which is followed by an op- 
tically Seyfert phase in massive galaxies with relatively 
high SFR in their host galaxies. 
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